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An excitation method for MEMS devices with planar electrodes is described. The stationary 
part of the device (the stator) consists of three electrode arrays arranged in the ABC ABC order. 
A, B, and C carry time-independent potentials and together form a spatially-periodic electrostatic 
profile. The moving part of the device (the translator) has two electrode arrays ababab, with a and b 

~f^ , carrying time-dependent out-of-phase voltages. When the frequency of the time-dependent voltage 

is close to the natural frequency of the spring-mass system, the translator is driven into resonance. 
By adjusting the spatial phase of the stationary profile, the driving force on the translator can be 

^\ maximized for any equilibrium position. Physical misalignment between the stator and translator 

resulting from imperfect fabrication can be corrected electrically. A dynamical equation describing 
translator motion is derived and analyzed for resonant and parametric driving. In both cases, the 
driving force depends on the translator displacement in a periodic fashion. Such nonlinearity of the 
driving force results in self-stabilization of forced oscillations. This property has implications for 

{T 1 } ■ the stability of vibratory gyroscopes. 
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I. INTRODUCTION 

In vibratory MEMS gyroscopes proof masses are driven into resonance to maximize the amplitude of a Coriolis 
force [l|. Amplitude or frequency variations of this motion directly translate into variations of the Coriolis signal. 
Gyroscopes with low angle random walk require stable primary oscillations insensitive to temperature, pressure and 
other environmental changes. One way to excite the proof mass is by using comb drives 043- ^ n this method, the 
driving force is independent of the proof mass position and the oscillation amplitude is determined by the balance 
between external, elastic and dissipation forces. Changes in the quality factor Q due to environmental variations 
directly affect the amplitude, which in turn causes the scale factor to drift. In practical devices, these effects are 
typically mitigated by an amplitude control loop |l|. 

The intrinsic Q-sensitivity of the resonator can be reduced in nonlinear drives if the driving force amplitude decreases 
with displacement. (Note that in common gap-closing drives the force increases with displacement.) An example of 
such a system would be the magnetic pendulum in crossed steady and oscillating magnetic fields [9| or the gravitational 
pendulum driven by a horizontal harmonic force (lCf . Since only the tangential component of the external force affects 
the motion, the driving force decreases away from the equilibrium position. As the resonator is excited, energy transfer 
slows down and approaches zero when a stationary amplitude is reached. As long as dissipation is not dominant, the 
final stationary amplitude will be determined by the restoring and external forces only, thus eliminating dependence 
on Q. 

In this paper, we describe another nonlinear driving method, referred to as three-phase driving, which has a similar 
property of self-stabilization. Unlike pendulum, the proof mass undergoes linear motion and is excited by planar 
electrodes deposited on two flat surfaces separated by a micron-scale gap. The electrode arrangement is similar to 



that of three-phase accelerometer [11| , but in this case electrostatic interaction between the stationary and moving 
electrodes are used to excite mechanical motion. The electrode geometry is defined in section UU capacitance matrix 
analyzed in section IIII1 and electrostatic force derived in section HVl The nonlinear equations of motion are derived 
and solved in section IVII for resonant excitation and in section IVIII for parametric excitation. 

Another important feature of three-phase driving is misalignment tolerance. By adjusting the voltages of stationary 
electrodes one can always tune the maximum of the driving force to coincide with the mechanical equilibrium of the 
proof mass. Any misalignment between the moving and stationary parts caused by imperfect fabrication can be 
compensated electrically. This property is derived in section [V] 

II. ELECTRODE GEOMETRY 

Both stationary and moving electrodes are periodic arrays of parallel strip lines extended along the y-axis. It is suffi- 
cient to consider only the two-dimensional (xz) cross section of the system, which is shown in figure [T] The stationary 
electrodes are split into three different groups A, B, and C, and arranged in a periodic sequence ABC ABC ABC . . . 
Physically, all the electrodes are identical, with equal widths and equal gaps in between. The center-to-center distance 
between two nearest electrodes of the same group is L, which is the spatial period (pitch) of the array. The smallest 
center-to-center AB distance (equal to the smallest BC and CA distances) is L/3. All electrodes that belong to the 
same group are electrically connected and carry the same time independent potential. The potentials are set according 
to the three-phase driving rule: 

<f> A = V coa (0), (1) 

4>b=Vcob(o+?£\, (2) 

<t>c = Vcos(e+^p\ =Vcos^- yY (3) 

The amplitude V and phase 9 are two adjustable parameters. In practical devices, V and 6 are set by external 
electronics. Collectively the A, B, and C electrodes will be referred to as the stator. 

The moving electrodes are split into two groups, to be denoted a and b, and arranged in the sequence abababab . . . 
Again, the electrodes are physically identical and equidistant. Their width and gap are in general different from those 
of A, B, and C, but the total spatial period is equal to L. The center-to-center distance of the nearest ab pair is 
L/2. Collectively, the ab electrode array will be referred to as the translator. In three-phase driving, a and b carry 
out-of-phase time-dependent voltages 

cj> a (t) = U cos (wt) = U(t), (4) 

4> b {t) = U cos (u>* + tt) = -U{t). (5) 
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FIG. 1. The basic electrode geometry. A, B, and C are stationary. The translator electrodes a and b can move along the x 
direction, x is referenced from the position in which the center of a aligns with the center of A. The electrodes are supported 
by dielectric substrates with a static dielectric constant of 12.1 (single-crystal silicon). The middle gap has a dielectric constant 
of 1.0. Notice the periodic boundary conditions along the a;-axis. Typical dimensions are (all in microns): array pitch L — 24.0, 
bottom dielectric thickness = 10.0; bottom electrodes width = 4.0; bottom electrodes height = 0.2; top dielectric thickness = 
10.0; top electrode width = 6.0; top electrode height = 0.2; the stator-translator gap = 1.6. 

The amplitude U and the angular frequency w are set by external electronics. The overall phase of U(t) is unimportant 
and will be set to zero. 

The rest of the device is represented by the mechanical support, which is electrically insulating, and by the frame 
electrode /. The latter is always electrically grounded, <f>f = 0. The frame electrode is assumed to be uniform along 
the x axis. 

To describe translator motion, a coordinate reference is needed. In this paper, x — corresponds to the center of 
an a electrode coinciding with the center of an A electrode, as indicated in figure [T] Assuming linear elastic force and 
linear dissipation, and neglecting static friction and stochastic noise, the equation of motion for a translator mass M 
reads 



M x + bx + kx — 



dW(x) 
dx 



(6) 



Here W{x) is the quasi-electrostatic energy of the system, whose x-derivative is the external force acting on the 
translator. An expression for W{x) under the three-phase rules is derived in the next two sections. 

III. CAPACITANCE MATRIX 

The device under study is a quasi-stationary system of conductors with defined electrostatic potentials. The full 
electrostatic energy is given by 



N N 



(7) 



1,3 =1 



where N — 6 is the total number of conductors. (N — A, B, C, a, b, /.) The symmetrical capacitance matrix dj is 
defined as 



N 

8=1 



'3 ■ 



(8) 



In three-phase driving, the capacitance coefficients are functions of the translator position x. As a result, the electrode 
charges are also position-dependent, qi = qi(x). 
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TABLE I. Symmetry properties of the capacitance functions. "Even" means dj(—x) — dj(x). "=" means this relationship is 
chosen as the definition of the generating function in question. For example, Gt{x) is defined as Co,a(x). 

In general, the energy (JT)) comprises 21 different contributions and depends on N(N + l)/2 = 21 capacitance 
functions. Since the frame electrode is assumed to be always grounded, <j)f — 0, the number of terms is reduced to 
15. In addition, translation, reflection, and permutation symmetries within the A, B, C, and a, b groups reduce that 
number to just 8 irreducible functions. The symmetry properties are summarized in table HI 

The capacitance coefficients Cy can be calculated numerically using finite-element software such as COMSOL 
Multiphysics. Following the definition ©, the potential of one conductor is set to 1 V and the others are grounded. 
Then a two-dimensional Laplace equation is solved and the charge induced on all the conductors is calculated. Notice 
that the solution domain can be reduced to one irreducible unit cell < x < L by imposing periodic boundary 
condition <p(L, z) — (p(0, z) for all z. The procedure is repeated for multiple translator shifts x to construct the entire 
capacitance functions. Several exemplary functions are shown in figure [2j The electrode array pitch is L = 24 /xm 
and the stator-translator gap is 1.6 /im. Other parameters are specified in the caption of figure Q] It can be observed 
that the functions shown in figure [2] satisfy the translation and permutation symmetries listed in table |U 

All eight generating functions are periodic and even. As such, they can be expanded in cosine Fourier series 



Gi{x) = yj G^cosfc„a; 

n=0,l,... 



(9) 



where fc„ = 27rn/L, and the prime at the sum sign means the n = term is taken with weight 1/2. All Fourier 
coefficients G„ are real numbers and have the same dimensionality as the original Gi(x). If the function itself is 
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FIG. 2. The capacitance functions Cij(x) for the array pitch L = 24 /jm and stator-translator gap of 1.6 /im. Other parameters 
are the same as in Fig. [T] The capacitances are per unit cell of the array, per unit length in the y direction (1 m), and given in 
units of £o = 8.854 • 10~ 12 F/m. Note different scales in different panels. The translation and permutation symmetries should 
be compared with table |T] 

u\ ... . . i I 

known then G„ can be found by numerical integration. The Fourier coefficients for the parameters of figure [U and 
functions of figure [2] are summarized in table [TTJ The coefficients decay rapidly with n. 

The dynamics of the translator is for the most part determined by the generating function Gt(x), see section [TV] 
This function is smooth and its Fourier expansion converges very rapidly with the number of harmonics included. 
Already the first three coefficients reconstruct the function with a better than 1% accuracy. Five coefficients are 
sufficient to accurately represent G*?(x) at all x. 



IV. ELECTROSTATIC ENERGY AND FORCE 

Referring to equation ([7]), the total number of energy terms is 21. However, since </>/ = 0, only 15 terms remain. 
Next, one makes use of the three-phase rules Q and ([5]) implying that the a and b potentials are always out-of-phase. 
Expanding the double sum, one obtains 



W=-U 2 (t)(C aa -2C ab + C bb ) 

+7j (Caa<Pa + Cbb4> 2 b + Ccc4> 2 c) + {Cab^a'Pb + Cac4>a4>c + Cbc4>b4>c) 

+ U(t) [MCaA - C bA ) + MCaB - C bB ) + 4>c{C aC - C bC )\ . 



(10) 



The first line here involves only the translator electrodes and therefore can be interpreted as a translator self-energy. 
The second line involves only the stationary electrodes and can be interpreted as a stator self-energy. The last line is 
a stator-translator interaction energy. Note that the translator and stator self-energies are in general functions of x 
and as such may contribute to the dynamics. 

The capacitance functions are now expressed via five generating functions according to the rules of table HI and the 
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TABLE II. First four non- vanishing Fourier coefficients for the eight generating functions. The values are given in units of eo. 
The geometry parameters are given in the caption of figure [T] The zero-order coefficients of the frame functions Gg , Gq and 
Gg are unphysically large because the dielectric thicknesses are unrealistically small (10 /an). These coefficients do not affect 
the electrostatic force on the translator. 



stator potentials are substituted from equations ([TJ-Q. The result is 



W(x) = W tT (x) + W st diag (x) + W st off -diag(a:) + W int (x) 



U 2 (t) [Gtix) + Gi(sc + 1/6) - 2G 2 {x)} 



+ -V 2 [G A {x) cos 2 (0) + G 4 (x - L/3) cos 2 (0 + 2tt/3) 

+ Gi(x + L/3) cos 2 (0 - 2tt/3)] 

+ V 2 [G 5 (x) cos (0 + 2tt/3) cos {9 - 2tt/3) 

+ G 5 (x - 1/6) cos (6>) cos (6> + 2tt/3) 
+ G 5 (a; + L/6) cos (9) cos (0 - 2tt/3)] 

+ VU(t) {cos (6») [G 7 (x) - G 7 (x + 1/2)] 

+ cos (0 + 2tt/3) [G 7 (x - L/3) - G 7 (.t + L/6)] 
+ cos (6» - 2tt/3) [G 7 (x + L/3) - G 7 (x - L/6)]} . 



(11) 



Next, Fourier expansions should be substituted in place of the capacitance functions and the resulting expressions 
simplified. The four contributions will now be considered separately. 

Translator self-energy. When summing the Fourier series for Gi(a;), odd terms cancel while even terms double. As a 
result, the translator self-energy assumes the form 



OO I 

W tT (x) = U 2 (t)J2[G ( <£-G, 



n=0 



cosfc6„a; . 



(12) 



Diagonal stator self-energy. The diagonal stator self-energy involves generating function G&{x). Upon substitution of 
the corresponding Fourier expansion, it is convenient to consider three sets of harmonics separately: (i) 0, 6, 12, . . . = 
6n; (ii) 2, 8, 14, ... = 2 + 6n, and (iii) 4, 10, 16, ... = 4 + 6n. In all three cases, a combination of trigonometric functions 



is simplified to a concise expression. The final result is 



W st diag (x) = - V 2 \2 V G$ cos k 6n x 

8 I n=0 

OO 

+ ^2 I G 2+6« C0S ( k 2+6nX + 29) + (?4+ 6 „ cos (fc 4+6 „a; - 29) 



(13) 



Off-diagonal stator self-energy. Again, three sets of G$(x) harmonics are treated separately. After some algebra, the 
off-diagonal stator self-energy is 
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n=0 
(5) 



2+6 „ cos (fc 2+6 na; + 20) + G\l 6n cos (fc4+6n£ - 20) 



(14) 



The structure of equation (Ti"4"|) is similar to that of the diagonal stator energy, equation (IT31 . 

Stator-translator interaction energy. In the last term of equation dill) harmonics with even n cancel while those with 

odd n double. The interaction energy becomes 

OO 

W int (a:) = 2VU{t) ^ G^ [cos k n x cos (0) 



n=l,3, 



+ cos (/c n a; — 27rn/3) cos (0 + 2ir/3) 
+ cos (k n x + 2irn/3) cos (9 — 2ir/3)) . 

Transformations are convenient to perform separately for the following groups of harmonics: 1, 7, 13, . 
3, 9, 15, . . . = 3 + 6n, and 5, 11, 17, . . . = 5 + 6n. The final result is 

OO 

W int (x) = 3VJ7(t) Y, { G i+6u cos {k 1+6n x + 9) + G^ 6 „ cos (k 5+6n x - 9)} . 



= 1 



(15) 
- 6n, 

(16) 



Note that harmonics with n = 3, 9, 15, . . . have vanished. This is a consequence of three-phase rules ([TJl-Q. 

Total electrostatic energy is given by the sum of equations (|12p. (|T3"j) . (|14p and l|16p. The electrostatic force can be 
found as a negative derivative of the energy with respect to translator position x. Performing differentiation one finds 
(recall that k n = 2im/L): 



F el _ 8t (x) = U 2 (t) £ k 6n [<» - G&> 
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OO 

+3VJ7(t) J2 { fc i+6«G^ 6n sin (fc 1+6 „x + 9) 

+ fc 5+6n G^ 6 „ sin (fc5+6«a; - 9) } ■ 



n=0 



(17) 



Formula (fTTf is general. Based on the results of finite-element modelling summarized in tableHH the general expression 
can be reduced to a handful of the most relevant terms. The following observations are made, (i) The translator 



self-force begins with coefficients Gg and Gg . According to table [TT1 both coefficients are of the order 0.01 or 
less, and hence can be neglected. The translator self-force does not contribute to its own dynamics, (ii) The stator 
self-force contains even harmonics of functions G±{x) and G${x). Only the n = 2 coefficients are larger than 0.04; 
the rest can be omitted. The stator self-force does contribute to the translator dynamics, (iii) The interaction force 
starts with coefficients G\ and G 5 . The latter is of the order 0.0003, and should be omitted. Thus although the 



8 

main driving function Gt(x) is highly anharmonic, its effect is represented by the sole coefficient C{ with a very high 
degree of accuracy. Putting all of this together, one arrives at a truncated electrostatic force in the following form: 

F d _ rt (aO = \v 2 [fa [g 2 4) + 2G 2 5) ] sin (fax + 29)} 

+ SVUcoautikiGP sin (fax + 6)\. (18) 

This expression can be used for all practical purposes. Including other small terms would have exceeded the accuracy 
of the numerical solution and of the model itself. The first term in the force originates from the stator self-energy. It 
is present even if translator electrodes are grounded. This force has no explicit time dependence and can be thought 
of as a correction to the elastic force acting on the translator. This correction is proportional to the square of the 
three-phase stator amplitude V 2 . The second term in equation (JT5J) originates from the interaction between the stator 
and the translator. It oscillates in time with a frequency ui and is the main driving force of the device. This force is 
nonzero only if both the stator and translator are at potentials prescribed by the three-phase rules. 

V. SELECTION OF 9 

Choosing optimal 6 is now discussed. In the presence of stator-translator misalignment, the translator is perma- 
nently shifted by a distance xq from the designed position. The elastic energy of a linear spring is 

-2 

Al" ft Hi 1* 

W e i ast ic(a;) = -(x-x ) = — , (19) 

where x — x — xq is the new independent variable. Since time derivatives of x are the same as of x, equation of motion 
([6]) becomes 

Mx + bx + nx = F c \- st (x + x ) . (20) 

Consider now the main driving force, which is the second term in equation (|18[) . Its x-dependence is defined by the 
factor 

sin (fax + 9) = sin (fax + kiXQ + 9) . (21) 

In forced oscillations, it is desirable that the external force is maximal at the equilibrium position x = where the 
velocity is maximal. This condition maximizes the mechanical power absorbed by the translator. Expression (1211) is 
maximized by choosing 9 such that 



fax + Q = \- (22) 



ZfaVUG ( p cos (fax) cos ut (23) 



Then the driving force 



is maximal at x = 0, as desired. The relation (|22[) expresses an important property of three-phase driving. By 
adjusting three phase 9 electronically one can always cancel the effects of physical misalignment, correct for process 
variations, and maximize the driving force. 

The stator-stator induced force is defined by the following factor [the first term in equation (|T5)) ; note that fa = 2 fa]: 

sin (fax + 29) = sin [fax + 2(kiXQ + 9)} = sin (fax + n) = — sin (fax) . (24) 

Collecting all the terms one obtains the final equation of motion 

Mx + bx + kx+ | k 2 V 2 \G { 2 ] + 2G 2 5) | sin (fax) = 3 faVUG { p cos (fax) coswi , (25) 

which is a central result of the paper. The last term on the left does not depend explicitly on time, and is interpreted 
as an electrostatic contribution to the translator restoring force. Note that it has the same symmetry as the elastic 
force: it is odd in x. At small vibration amplitudes the electrostatic force simply renormalizes spring stiffness k. At 
large amplitudes, it brings about new effects. Both the electrostatic force and the driving force are periodic functions 
of the translator displacement x. Thus the equation of motion (|25p is highly nonlinear. 



Parameter 


e 


5 


g 


n 


Value 


0.00076 a 
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0.036 c 


«l d or «2 e 



TABLE III. Typical values of the dimensionless parameters (I27ll - (|30| l for the geometry of figure [T] The estimates assume 
the translator wafer thickness 125 fxm, silicon mass density 2330 kg/m 3 , natural frequency ujo = 2-k ■ 6500 1/s, stator voltage 
amplitude V = 8 volts, and translator carrier amplitude U — 12 volts. a Estimated from ringdown data. Q — l/(2e) = 658. 
b Estimated from equation (I28[) and table [II] c Estimated from equation ()29[) and table [Tl] d Resonant driving. c Parametric 
driving. 



To facilitate subsequent analysis, it is convenient to rewrite the equation of motion in scaled variables, (i) Introduce 
dimensionless displacement £ and time r: 



£ = k\x = 



2-kx 



: UJot = t 




(26) 



where loq = ^/k/M is the translator's natural frequency at small amplitudes. Hereafter, the derivative with respect to 
r will be denoted by a prime rather than a dot. (ii) Divide the equation of motion by (Mu>qL)/(2tt). (hi) Introduce 
new dimensionless parameters 



2Mw ' 



(27) 
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(28) 



ZVUG ( P 

Mcu 2 



(29) 



The equation of motion (1231) becomes 



Cl 






i" + 2ei' + £ + S sin (2£) = g cos (£) cos fir. 



(30) 



(31) 



The four parameters have the following meaning, (i) e = 1/(2(5) determines the rate of energy dissipation in the 
system; it is an inverse of the quality factor Q. (ii) 5 is the measure of electrostatic force; it is proportional to V 2 ; 
(iii) g is the dimensionless amplitude of the main driving force; (iv) fl is the dimensionless external frequency that 
characterizes frequency mismatch between the driving force and the natural frequency. Typical parameter values are 
given in table [ 



VI. RESONANT DRIVING 



The most interesting feature of the equation of motion ([3~T|) is the nonlinear driving force g cos (£) cos fit. Nonlinearity 
of such type have been studied in relation to chaotic behaviours of driven magnetic pendulum 9] , horizontally driven 
gravitational pendulum [id ], elliptically driven gravitational pendulum [I2l4l4| . and charged particles in the field of 
plane waves [15|-[l7|. In those physical systems, the functional form of the restoring force is linked to the functional 
form of drive nonlinearity. For example, the restoring force of the gravitational pendulum is proportional to sin£ and 
has the same origin as the cos£ of the driving force. In the system under study, the electrostatic correction (5sin2£ is 
similarly related to the driving force. However, the other part of the restoring force, £, comes from elastic flexures, 
i.e., from a different part of the system. In fact, other terms can be added to the restoring force, for example the 
Duffing term £ 3 , without any need to modify the driving force. The independence of restoring and driving forces adds 
to the richness of this dynamical system. 
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The primary purpose of the present work is to study the nonlinear dynamics at small g when oscillations are stable 
and no chaotic motion is observed. In addition, the electrostatic force will be neglected in the following analysis, 
despite its definite presence in practical devices. Setting 5 = 0, the equation of motion (|3T]) reduces to 



£" + 2e£' + £ = g cos (£) cos Or . 



(32) 



In this section, resonant driving with !1«1 will be considered. The goal is to derive resonance curves by applying 
the method of slowly changing amplitudes. A solution to equation (|3"2")l is sought in the following form 

£(t) = cr(r) cos [Vlr + ip( T )] . (33) 

In the stationary regime, the amplitude a and phase ip are time- independent. Substitution in equation (|32j) yields 

(l - Vl 2 ) a cos {Vlr + ip) - 2e Via sin (Vlr + tp) = 



— g cos (Vlr) cos {a cos (Vlr + ■;/>)}. 

Shifting the time origin by tp/Vl, the equation is brought to the form 

(l - Vl 2 ) a cos {fir) - 2e Via sin (Vlr) = 
= g cos (Qt) cos "0 cos {a cos (fir) } + g sin (Or) sin -0 cos {a cos (fir) } 



(34) 



(35) 



By multiplying with cos (fir) and sui($7t) and averaging over the oscillation period 2n/Vl, one obtains a pair of 
equations 



(1 - Vl 2 ) a = , 9 cos ip{J (a) - J 2 (a)} 
—2eVla — gsinip{J (a) + J2(c)} 



(36) 
(37) 



where J n (a) is Besscl function of order n. Applying the identities Jo (a) — J<i (a) = 2 J[(cr) and Jo(c)+ J2(cr) = 2Ji(<t)/<t, 
equations (f36 |) - (f37|) are rewritten as 



1-^ 



-2eO 



2J{(a) 

„2 



L2Ji(a) 
Squaring and summing these equations eliminates phase tp: 



gcosip 
g sin ip . 



'l-Vt 2 



2c 



\aJ[(*Y 


n| = 





It is convenient to introduce renormalized force and dissipation parameters: 

g r =g-[2Jl(a)}, 
-aJ[(a) 



,h{a) 



(38) 
(39) 

(40) 

(41) 
(42) 



The subscript r here stands for "resonance" . The most important features of g r and e r is that they tend to zero 
as the derivative of the n = 1 Besscl function approaches its first zero or, equivalently, as Ji(cr) approaches its first 
maximum. This occurs at a — 1.85. At small amplitudes, J\{o <1)k c/2, J[(<t <1) w 1/2, and g r — > g, e r — > e. 

Equation (|40|) establishes a relation between the drive frequency Vt and the amplitude a. It is easier to express Vt 
vs. a rather than vice versa: 



n r , 



cs. drive 



f/r 



(l-2e?)± 1/^-43+43 



(43) 



This resonance function is shown in figure [3] as the thick and thin solid lines. The thick line corresponds to weak 
damping e = 0.00076. At this value, the amplitude of forced oscillations increases until the renormalized force g r 
starts to decrease. The two processes balance each other and the stationary amplitude converges to 1.85. It must be 
emphasized that the amplitude would have been much larger than 1.85 in the usual case of linear driving at the same 
damping. Such a curve is shown in figure [3] by the dashed line. 
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e = 0.00076, linear drive 

e = 0.00076, nonlinear drive 

£ = 0.02, nonlinear drive 




1 1.1 

Frequency Q. 



FIG. 3. The resonance curve (|43[) corresponding to nonlinear drive and linear spring (thick and thin solid lines). The parameters 
are g = 0.036, 8 — 0, e = 0.00076 (the thick line) and e = 0.02 (the thin line). The dashed line is the linear drive resonance 



curve for e = 0.00076 and g = 0.036. 



One concludes that at weak damping the amplitude of forced oscillations is no longer limited by damping as in 
conventional linear resonance. Instead, it converges to a constant value of 1.85. In dimensional units, see equation (|26[) . 
it corresponds to 



1.85 

~2^T 



L«0.294L. 



(44) 



Thus three-phase driving possesses an important property of self-limitation. The stationary amplitude is weakly 
dependent on the force magnitude g and is directly proportional to the spatial period of ABC array. Another 
property of the nonlinear force is that the stationary amplitude is weakly dependent on damping rate e as long as the 
latter is below a critical value. 

At larger e (small Q), damping becomes so strong that at some point it starts to limit the amplitude. It happens 
when the tip of the resonance curve detaches from a = 1.85. In other words, the inner root in equation ((43]) becomes 
zero at a — 1.85. Expanding this condition and neglecting if relative to e£, one obtains 



9- 



J,(a) 



0.169g. 



(45) 



<r=1.85 



For g = 0.036, this yields e cr = 0.0061 (and the critical quality factor Q cr = 82). This is about 8 times larger than 
e = 0.00076. Thus the parameter set of table Hill is well within the pitch-limiting regime. An example of damping- 
limiting regime, e = 0.02, is shown in figure[3]by the thin solid line. In this case, the stationary amplitude only grows 
to 0.82 at resonance, a value well below the pitch limit of 1.85. 

The above analysis is supported by direct solutions of the equation of motion (|32p . Figures U and [5] compare time 
evolution of the translator position in the pitch-limited and damping-limited regimes. In figure [4j damping is weak, 
e = 0.00076, and motion is underdamped. Initially, the amplitude grows fast and quickly outgrows the stationary 
value of 1.85. Then the driving force effectively changes sign and becomes a stopping force. As a result, the amplitude 
drops below 1.85. Such swing cycles continue for some time until their intensity subsides and the amplitude settles 
at the stationary value. In contrast, in figure [5] e = 0.02, and the motion is overdamped. No amplitude swings are 
observed. Instead, the amplitude rises monotonically and smoothly approaches a stationary value of 0.82. 



VII. PARAMETRIC DRIVING 



In resonant driving considered in the preceding section, the external force was set in phase with translator velocity x 
which led to the phase condition (|22[) and the equation of nrotion (1311) . This condition maximizes the power transferred 
to the oscillator. In parametric driving, the external force is set in phase with translator displacement x. That is, 
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2000 4000 6000 8000 10000 

Dimensionless time x 



FIG. 4. Numerical solution of equation (|32|) for Q, — 1.0, e = 0.00076 and g — 0.036. After swings oscillations, the amplitude 
settles at 1.85. The initial conditions are £(0) = £'(0) = 0. 




100 150 200 250 300 

Dimensionless time x 



400 



FIG. 5. Numerical solution of equation (|32[) for Q = 1.0, e = 0.02 and g — 0.036. The amplitude gradually approaches the 
stationary value of 0.82 with no swings. The initial conditions are £ (0) = £'(0) = 0. 



the main driving force is zero when x = 0. From this and equation (|21[) follows a different, parametric three-phase 
condition 



The main driving force is [the last term in equation (j!8[) ]: 

(7) 

± 3kiVUG\ sin(fcix) cosujt 



(46) 



(47) 



The overall sign is plus or minus depending on whether zero or n is chosen. As a matter of fact, the choice is irrelevant 
since the sign can be flipped by shifting the time origin by tt/oj. It is convenient to choose the negative sign. The V 2 
term in the truncated force (fT5|) is defined by the factor 



sin (k2X + 28) = sin \k^x + 2(kiXo + 9)] = sin {k^x + {0, (27r)}) = sin (k2x) 



(48) 



Substituting the truncated force into the full equation of motion and switching to dimensionless units as in section IVT1 
one obtains 



i" + 2e£' + £ - 6 sin (2£) + g sin (£) cos fir = 0. 



(49) 
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Comparing to resonant driving, this equation of motion has a different sign of the S term and a sine instead of cosine 
nonlinearity in the driving force. The sin£cosili nonlinearity has been extensively studied over decades in relation to 
parametrically excited pendulum }18l - t3ll [ . Again, the system under study is more general. Its elastic force is decoupled 
from the driving force and as such can be of arbitrary form allowed by the physics of elastic springs. 

The method of slowly changing amplitudes can also be applied to parametric driving. Limiting again consideration 
to S = 0, the equation of motion (|49|) reduces to 



£" + 2e£' + £ + g sin (£) cos fir = 
The solution is sought in a single-harmonic form 

£(t) = ct(t) cos [wt + iP(t)] , 



0. 



(50) 



(51) 



where ui f» 1 is close to the natural frequency of the proof mass. The external frequency is set to 51 = 2to. In the limit 
of stationary oscillations, a are ip are treated as constants. Substituting in equation (|50p and shifting the time origin 
by ip/u yields 



(l — to J a cos lot — 2euicrsmujT — 
g cos (2tp) sin {a cos wt} cos (2ut) + g sin (2tp) sin {a cos wt} sin (2wr) 



(52) 



Multiplying by cos (lot) and averaging over time, then multiplying by sin (wt) and averaging, and finally applying the 
identities Ji(cr) — Jz(cr) = 2J' 2 (a) and Ji(er) + J3(er) = 4J 2 (e>-)/o', one obtains 



( t(1- W 2 ) = [2J^( ( t)]« ? cos(2^), 
gsm (zi/y). 



— 2ea;cr = 

By introducing new renormalized force and dissipation parameters 

.g p = .g-[2^(a)]=. g .[J 1 (a)-J 3 (a)] 



aJ 2 (a) 



2J 2 (a) 



J, (a) - J 3 (a) 
Ji(a) + J 3 (a) ' 



(53) 
(54) 

(55) 
(56) 



where subscript p stands for "parametric" , the equations are cast in the form 

a (l -co 2 ) =.9 p cos(2V0 , 
—2e p toa = g p sin (2-0) . 



(57) 
(58) 

This is the same functional form as equations (|38[) and (|39[) but with g and e replaced with g p and e p . Accordingly, 
the resonance curve is given by the function similar to equation (|43|) : 



^param. drive 



2e 



^2 



^/|-^ 



•4g*. 



(59) 



Both g p and e p tend to zero at the first maximum of the n = 2 Bessel function, which occurs at a = 3.05. By analogy 
with the resonant case, one expects that 3.05 will be the maximal possible stationary amplitude of parametric drive. 
In real units 



3.05 
~2^~ 



0.485 L. 



(60) 



One concludes that parametric driving results in about 65% larger amplitude than resonant driving (0.485 L vs. 0.294 

Another important feature of g p is that it tends to zero at small a <^ 1 as ex (c/2). Accordingly, the ratio g p j o 
tends to a constant at small amplitudes. As a result, stationary oscillations exist only within a small frequency 
interval around the natural frequency to = 1. This is in sharp contrast with resonant driving. The resonance curve 
(f59"| is plotted in figure [5] Its shape is distinctly different from that of resonant driving. Parametric oscillations exist 
within a small frequency window. The self-limitation property of parametric drive is confirmed by direct solution of 
the equation of motion (|49[) . Figure [7] shows time evolution of the translator amplitude in the underdamped regime, 
e = 0.00076. The amplitude exhibits slowly varying swings, similar to figure 21 before settling at a stationary value 
of 3.05, in accordance with equation ([6"tT|) . 
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FIG. 6. Parametric resonance curve (|59[) for e = 0.00076 and g = 0.036. External driving frequency is Q = 2uj. 




2000 4000 6000 8000 10000 

Dimensionless time x 



FIG. 7. Numerical solution of equation (|49[l for e 
£'(0) = 0.1. The stationary amplitude is 3.05. 



0.00076, 5 = 0, g = 0.036, and fi = 2.0. Initial conditions are £(0) = 0, 



VIII. 



SUMMARY 



In this paper, an excitation method for MEMS devices with planar electrodes has been presented. The key ingredient 
is the periodic electrostatic profile generated by electrode array ABC ABC . By independently adjusting voltages at 
A, B, and C, one can arbitrarily move the potential profile along the x-axis thereby compensating any fabrication 
misalignment between the stationary and moving parts of the device. A time dependent voltage applied to the 
translator array ababab interacts with the stationary profile and either drives the translator into resonance or induces 
parametric excitation. 

Using planar electrodes to induce sliding motion necessarily involves fringe fields. The underlying electrostatic 
problem is complex and requires careful analysis. Various symmetries of the capacitance matrix have been analyzed 
and summarized in table HI The capacitance coefficients have been computed by a finite-clement numerical method 
as functions of the translator position for a realistic electrode geometry. The results are shown in figure [5J Fourier 
coefficients are given in table [TTJ General expressions for quasi-electrostatic energy and force have been derived in 
sections IIIII and IIV1 respectively. Based on the analysis of Fourier coefficients, a truncated electrostatic force (|18p 
has been derived. In addition to the expected main driving term, the force contains a stator-induced self-force that 
affects the translator dynamics. The basic equation of translator motion has been derived for both resonant driving, 



15 

equation (J3TJ) . and parametric driving, equation (|49|) . 

Spatial periodicity of the electrostatic potential results in a periodic dependence of the driving force on the translator 
displacement. It leads to self-stabilization of forced oscillations. For weak dissipation, the amplitude of forced 
oscillation stabilizes at 0.294 of the electrode array pitch in the case of resonant driving, and at 0.485 of the array 
pitch in the case of parametric driving. Self-stabilization has been confirmed by direct solution of the equations of 
motion in time domain. 

Self-stabilization can be useful in applications where amplitude stability is required, for instance, in high- 
performance vibratory MEMS gyroscopes. Temperature, pressure and other environmental variations do not affect 
the array pitch to the same degree as they affect quality factor Q. Pitch limited nonlinear excitation methods may 
provide a higher degree of stability than dissipation limited ones. 
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